Steady-state crystallization of Rydberg excitations in an optically driven lattice gas 
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We study resonant optical excitations of atoms in a one-dimensional lattice to the Rydberg states 
interacting via the van der Waals potential which suppresses simultaneous excitation of neighboring 
atoms. Considering two- and three- level excitation schemes, we analyze the dynamics and stationary 
state of the continuously-driven, dissipative many-body system employing time-dependent density- 
matrix renormalization group (t-DMRG) simulations. We show that two-level atoms can exhibit only 
nearest neighbor correlations, while three-level atoms under dark-state resonant driving can develop 
finite-range crystalline order of Rydberg excitations. We present an approximate rate equation 
model whose analytic solution yields qualitative understanding of the numerical results. 
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Strong, long-range interaction between Rydberg atoms 
[l| have positioned them as promising systems for quan- 
tum information processing which motivated con- 
siderable experimental progress in preparing and study- 
ing such systems [5l-ll2| . Rydberg atoms are also interest- 
ing in the context of many-body physics: It was predicted 
that the long-range interaction leads to spontaneous sym- 
metry breaking and crystalline order in a continuous gas 
Hi-ill, formation of three-dimensional super-solids [lq . 
as well as to vortex lattices [ItJ or fractional quantum 
Hall states [HI in two-dimensional (2D) systems with 
artificial magnetic fields. Coupling Rydberg atoms to 
light can give rise to novel photonic states with highly 
non-classical correlations fl9l [20| | and unusual nonlinear 
spectroscopic features |2lL |22|. 

For regular arrays of coherently driven atoms, depend- 
ing on the strength and detuning of driving lasers, dif- 
ferent ground-state phases with crystalline order emerge 
[H-Iil. Rydberg-dressed atoms [ll, , i.e. ground- 

state atoms with a small admixture of Rydberg states, 
can also form crystalline structures with fractional fill- 
ings [29]. In the optically driven lattice gas [30|, where 
the number of Rydberg excitations is not conserved, 
the preparation of the ground state of the system re- 
quires careful consideration and an adiabatic preparation 
scheme has been proposed and analyzed in [3l|. A more 
natural approach to the formation of crystalline order of 
Rydberg excitations is to utilize the stationary state of 
a dissipative many-body system, which results from the 
interplay between continuous optical driving and sponta- 
neous decay. Moreover, the steady state is an attractor 
of the system's dynamics and is therefore stable against 
small perturbations. 

Here we study resonant optical excitations of Ryd- 
berg states of atoms using two- and three-level driving 
schemes. The van der Waals (vdW) interaction between 
the atoms leads to a Rydberg level shift and thereby 
blocks the excitation of an atom which is sufficiently close 



to an already excited Rydberg atom. Specifically, we 
study a one-dimensional (ID) lattice with strong nearest 
neighbor interaction between the atoms, employing nu- 
merically exact t-DMRG simulations. We show that two- 
level driving leads to at most short-range correlations of 
Rydberg excitation probabilities of atoms at neighbor- 
ing lattice sites. In contrast, for three- level atoms un- 
der the dark-state resonant driving, we find longer-range 
correlations and quasi-crystallization extending over sev- 
eral lattice periods. We derive an effective rate equations 
model, which can be solved exactly for the steady state 
yielding an expression for the correlation length. Finally, 
we estimate the influence of the tails of the vdW poten- 
tial inducing small level shifts experienced by the next to 
nearest neighbor atoms, which give rise to an upper limit 
of the correlation length. 

We thus consider a chain of N atoms trapped in a ID 
lattice potential of period a [Fig. DJa)] and examine one- 
and two-photon resonant excitations of atoms from the 
ground state \g) to the Rydberg state |r) [Fig. DJb)]. In 
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FIG. 1. (Color online), (a) Schematics of optically driven 
atoms in a lattice, (b) Starting from the ground state \g), 
the atoms are resonantly excited to the Rydberg state \r) 
either directly (left) or via resonant intermediate state |e) 
(right), (c) Illustration of the rate equations model with Ry- 
dberg blockade of neighboring atoms. 
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the two- level scheme, the transition \g) — » |r) is driven 
by a laser field of (effective) Rabi frequency Q gr . In 
the three- level scheme, the Rydberg state |r) is popu- 
lated from the ground state \g) via resonant intermedi- 
ate state \e) in the coherent population trapping (CPT) 
or dark-state resonance configuration [32[ with Rabi fre- 
quencies Q ge > Q er . The corresponding atom-field in- 
teraction Hamiltonians are given, respectively, by V J 2 — 



the master equation 



-ft(JV^ + H.c.) and V{ = -h(Q ge ai g ^Q er a^ re ^R.c.), 
where & = \p)jj{v\ are the transition operators for 
atom j. 

Spontaneous decay is described by Liouvillian terms in 
the equation of motion for the density operator p: Op = 

\{2UpL^ - {UW, p}), where V are the Lindblad 
generators. For the two- level scheme, 
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with being the (population) decay rate of Rydberg 
state |r). Although T rg is typically small (~ 10 4 Hz), it 
may be comparable to the Rabi frequency Q gr associated 
with either direct one-photon (UV) transition \g) — >• |r) 
with small dipole matrix element, or two-photon tran- 
sition via far off-resonant intermediate states. For the 
three- level CPT scheme, we take into account only the 
large spontaneous decay rate T eg (~ 10 T Hz) of the in- 



termediate exited state \e) via L J 3 
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rate T re of Rydberg state can safely be neglected in com- 
parison with Q er . 

For a single (isolated) two- or three-level atom under 
continuous driving, the steady-state population of the 
Rydberg state is given, respectively, by 



(a rr ) 



(<T rr ) 



2|fVI 2 +7, 2 9 + A 2 ' 



\n ge \ 2 (\n ge \ 2 + \n er \ 2 ) 



|^e| 2 + |^| 2 ) 2 + (7e 2 <,+2|^e| 2 )A2 



(la) 
, (lb) 



where 7^ = |r ^ v and A is the one- or two-photon de- 
tuning of the laser frequency from the \g) \r) transi- 
tion resonance. It follows from Eqs. (pQ) that the Rydberg 
state population (a rr ) of an atom is a Lorentzian function 



of A, with the width w — A /2|£^ r | 2 + 7^ for the direct 
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the CPT excitation schemes. 

Finally, pairs of atoms i and j interact with each other 
via the vdW potential [33[ V*^ w = ha l rr ^-a^ r , where 

dij = a\i — j\ is the interatomic distance. Once an atom 
i is excited to the Rydberg state |r), it shifts the atom 
j out of the resonance by A = U/\i — j| 6 , where U = 
Cq/cl 6 . The excitation of atom j to state |r) is then 
blocked if A > w, which determines the blockade distance 
<i b ~ yCo/w. Throughout this paper, we assume that 
a < d\) < 2a, so that there is a Rydberg blockade only 
between the neighboring atoms. 

The density matrix p of the system of N atoms obeys 



[Hup] +Cip, 



(2) 



with the Hamiltonian %i = J2j ^1 1 A^K j 
Liouvillian L\p — J2j for / = 2 or 3 level atoms. 

Since the interaction of an atom with its immediate 
neighbor is much stronger than that with the atoms fur- 
ther apart, we truncate the vdW potential V*^ w to near- 
est neighbor (NN), V^ N = fia l rr U a J rr for i = j — 1 and 
V^n = otherwise. We have performed exact numer- 
ical integrations of Eq. (|2]) for small systems of several 
(N < 7) atoms interacting via the V v aw and Vnn poten- 
tials, and verified that they yield similar results for both 
two- and three- level excitation schemes. We will there- 
fore employ from now on the Vnn potential, returning 
to the corrections due to longer-range interactions of the 
Vvdw potential later on. 

We obtain the time evolution and the steady state of 
the full many-body density matrix Eq. (|2]) employing the 
t-DMRG method [3, HHj]. Our implementation follows 
the original proposal by Vidal [36]] , generalized to open 
quantum systems [33,[38|. Simulations for up to N ~ 10 2 
atoms are possible as the relaxation keeps the entangle- 
ment inside the system small. As a consequence, the 
usual limitations concerning the propagation time do not 
apply and the t-DMRG integration can be performed for, 
in principle, arbitrarily long times. Furthermore, since 
the stationary state is an attractor of the dynamics of 
the system, accumulated errors self-correct. We verified 
that for small systems the exact and t-DMRG solutions 
with the Vnn potential are indistinguishable. 

Results of the t-DMRG simulations for two- and three- 
level atoms in a realistically large lattice of length 30a 
with open boundary conditions are shown in Fig. [2j 
Since interactions between the atoms suppress the Ry- 
dberg excitations, the atoms at the boundaries j = 
and j = N — 1 having only one neighbor acquire the 
largest population of the Rydberg state |r). We also 
observe period-2 spatial oscillations of excitation proba- 
bilities (cr J rr ). For two- level atoms, the amplitude of these 
oscillations, as well as the average Rydberg state popu- 
lation is low, both for weak Q gr < T rg and moderately 
strong Q gr ~ T rg driving, Fig. [2ja) and (b). This is due 
to the fact that even for a noninteracting two-level atom 
under strong driving Q gr ^> T rg the excited state popu- 
lation saturates to (a rr ) —> \ [cf. Eq. ([Ta]) ]. In contrast, 
for a (noninteracting) three-level atom under the CPT 
excitation scheme, the population of Rydberg state can 
be very large, (a rr ) — >• 1, when Q ge ^> Q er [cf. Eq. (|lb|) ]. 
For a chain of strongly interacting atoms, we then observe 
a large amplitude of spatial oscillations of the Rydberg 
excitation probabilities (<r£ r ) extending over many lattice 
sites, Fig. Efc). 

Atoms excited to Rydberg states are typically no 
longer trapped in an optical lattice, although there has 
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FIG. 3. (Color online). Dynamics of Rydberg excitation prob- 
abilities for simultaneous application of driving lasers at time 
t = to all the atoms (thin red lines), and for sequential 
application — "sweep" — of the lasers to atoms j = 0, 1, 2, . . . 
(thick black lines). All parameters are as in Fig. [2jc) and 
time is in units of ~ — 7? — . 



FIG. 2. (Color online). Steady-state Rydberg excitation 
probabilities (a J rr ) of two- and three- level atoms in a lattice of 
length 30a (N = 31 atoms). Black dots (connected by dashed 
lines) show the t-DMRG solutions of Eq. ([2]) with the trun- 
cated potential V^ N 5 an d blue stars are the REs solutions for 
the hard-core (infinite NN) potential. For the two-level atoms, 
the parameters are U = 2Q gr , Y rg = Q gr (a) and Y rg = \Qgr 
(b). For the three- level atoms, U 

^er — 5^fi ,e (c) • 
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been some progress in trapping Rydberg atoms using 
ponderomotive potentials [39]. An important question 
is therefore, whether the steady-state crystalline order of 
Rydberg excitations can be reached during times short 
enough for the center-of-mass motion of the excited 
atoms to be negligible. In Fig. [3] we show the time evolu- 
tion of Rydberg state populations of atoms initially in the 
ground state \g). Note that applying the driving lasers 
to all the atoms at once yields the global stationary state 
of the system after times two to three orders of magni- 
tude longer than the steady-state equilibration time of an 
isolated atom. The reason for this is the initial formation 
of finite-size anti-ferromagnetic domains with dislocation 
defects in between, which require excessive healing times. 
But the steady-state is independent of the initial condi- 
tions and details of preparation, therefore faster schemes 
for preparing the global stationary state of the system 
may exist. Similar to classical crystals, it is indeed much 
faster to "grow" the Rydberg quasi-crystal by first apply- 
ing the lasers to the atoms at the lattice boundary and 
then gradually extending the irradiated region until the 
full Hamiltonian H3 is realized. In Fig. [3] we verify this 
intuitive approach. Since the crystalline order is rooted 
in the Rydberg blockade of the NN sites, a "sweep" ve- 
locity as fast as one lattice period a per single-atom equi- 
libration time r eg /(Q ge Q er ) can be applied to reach the 



steady-state configuration. 

We now describe an effective rate equations (REs) 
model, which yields an analytic solution for the station- 
ary state of the system. We are interested in the Rydberg 
excitation probabilities (<3^ r ) of atoms and their steady- 
state correlations (<5^ r <7^ r ), which allows us to disregard 
the coherences between the atoms. We restrict the in- 
teratomic interactions to the complete blockade of Ry- 
dberg excitations of NNs, assuming the potential V^j N 
with U —> 00. At each lattice site, we then have two 
incoherent processes: a state dependent pump with rate 
P and a deexcitation with rate D, 
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The ratios of the two rates k\ 

I = 2 and 3 level atoms are obtained from Eqs. (pQ) as 
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In this model, the density matrix p(t) obeys the equation 
of motion 

p = J2 (?LiP^ ~ {LpLi, p} + 2i\pL* - {L*Ll p}) . 

3 

(4) 

After sufficient relaxation time, the density ma- 
trix attains an essentially classical form p(t) = 
E {nj }MW)W where p({nj}) is the 

probability of configuration {rij} G (0, 1)^. Classical 
state-space dimension grows exponentially with TV, pre- 
cluding numeric integration for large systems. However, 
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the steady state of REs fulfills the detailed balance rela- 
tion 

P({ m j}) _ Ytjimj-nj) 



p({nj}) 



(5) 



States with the same number of excitations have equal 
weight and the partition function is given by Zn = 
El ^(M,iV)^, where Q(M,N) is the number of 
possible arrangements of M excitations on a lattice of TV 
sites. In ID, we have the analytic expression Q(M,N) — 
/iV-M+i^ wn j cn enables efficient calculation of all the 
steady state probabilities p({nj}). 

In Fig. [2] we compare the solutions of the REs model 
with the results of t-DMRG simulations, observing rea- 
sonable agreement, especially for the three- level excita- 
tion scheme of the atoms. In all cases, however, the decay 
of correlations is correctly captured by the REs model. 

For a large lattice N ^> 1, the partition function con- 
verges to Zn = cosh ^(2 + A^)^") > w ^ n which the Ry- 
dberg excitation probabilities at odd sites j (the bound- 
ary being at j = 0) are given by (<3^ r ) = e _1 ^ 2 ^^^(l + 
e -j/ v^)/2. The decay of spacial oscillations of probabil- 
ities {cr 3 rr ) is characterized by the correlation length 



(6) 



Note that in a translationally invariant system (i.e. in 
the bulk), the excitation probability is uniform in space, 
but the onset of crystallization is revealed by the density- 
density correlations, which can also be calculated in the 
REs model. We find that is indeed the correspond- 
ing correlation length. 

For the two- level excitation scheme of atoms, < 1, 
we then obtain the correlation length £ < a, i.e. shorter 
than the lattice spacing. Using a mean-field approxima- 
tion for the same system in the stationary state, Lee et 
al. [Ill predicted a phase transition to anti-ferromagnetic 
order (period-2 density wave). The reason for this dis- 
crepancy is that the mean field approximation is inade- 
quate in ID. We note that even in a 2D square lattice the 
phase transition to the antiferromagnetic state cannot be 
reached with the two- level atoms; the 2D extension of the 
REs approach leads to the classical hard-square model, 
for which a phase transition to the Neel ordered (checker- 
board) phase occurs at the critical value of ft cr i t ~ 3.7962 

For the CPT excitation of three- level atoms, one can 
tune ^3 by the Rabi frequencies Q ge and ^ e r, attaining 
large but finite correlation lengths £ > a. Further in- 
crease of £ is constrained by the validity of the NN inter- 
action approximation. We therefore consider the correc- 
tions originating from the hitherto neglected interactions 
of the atom with the next nearest neighbors (NNNs). To 
estimate the maximal effect of the NNN interaction, let 
us assume a strong driving, resulting in Rydberg exci- 
tations at every other lattice site (half filling). A non- 
blocked atom is then likely to be detuned by A' = A/2 6 



due to the interaction with the NNN Rydberg atoms. 
Expressing the NN detuning A = j3w in terms of the 
excitation linewidth w and using Eqs. (fib)) , we obtain a 
modified 
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To ensure almost perfect NN blockade, it is reasonable 
to take j3 w 10 and ~ 1% excitation probability for the 
blocked atom. In ID, the correlation length is there- 
fore limited to values of £ < 7a, due to the softness of 
the vdW interaction potential and the long wings of the 
Lorentzian excitation profile. In 2D, the phase transi- 
tion to Neel order may occur for three- level atoms, but 
since the effects of the NNN interactions will be more 
pronounced, the existence of a phase transition in 2D 
remains an open question. We finally note that for the 
dipole-dipole interaction V^ D ex Eq. ([7]) reduces to 

^3 < (8//^) 2 leading to a short correlation length £ < a. 

To conclude, we have found that the collective steady- 
state of an ensemble of three- level atoms in a ID lattice 
can exhibit quasi-crystallization of Rydberg excitations 
with a correlation length extending over many lattice pe- 
riods. In contrast, for two- level atoms even under strong 
driving, the correlations are only between the neighbor- 
ing atoms and the phase transition in 2D to Neel order 
cannot occur. Using the time-dependent DMRG simula- 
tions for several tens of atoms, we have shown that for 
uniform optical driving of all the atoms, long-range crys- 
talline order of Rydberg excitations is attained only after 
very long times. A sequential excitation of neighboring 
atoms by dynamically "sweeping" the lattice with the 
driving laser, can result in the same steady-state of the 
system in a much shorter time given by the single-atom 
equilibration time multiplied by the number of sites. We 
derived an effective rate equations model whose exact 
steady-state solution, being in good agreement with the 
numerical calculations, yielded explicit analytic expres- 
sions for the Rydberg excitation probabilities and corre- 
lation length. 
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